Embedding a chaotic signature in a periodic train: can periodic signals be chaotic? 
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We show how a chaotic system can be locked to emit a periodic waveform belonging to its chaotic 
attractor. We numerically demonstrate our idea in a system composed of a semiconductor laser 
driven to chaos by optical feedback from a short external cavity. The clue is the injection of an 
appropriate periodic signal that modulates the phase and amplitude of the intra-cavity radiation, a 
chaotic analogy of conventional mode-locking. The result is a time process that manifests a chaotic 
signature embedded in a long-scale periodic train. 
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Chaos is an important behavior of dynamical systems 
that arises in dire distinction to standard Newton-like 
behavior in a great variety of physical arenas [l] . Chaos 
has been observed, for instance, in electrical circuit, fluid 
dynamics and laser. Chaotic behavior is a phase of the 
dynamics, and it is an accepted notion that a system 
that is chaotic is not periodic, and viceversa. Periodicity 
and chaos are, however, closely related. For example, 
a stable system can be routed to chaos passing through 
a limit cycle, which represents a periodic trajectory in 
phase space. The regularity of a limit cycle is however 
a distinct signature that makes it easily distinguishable 
from a chaotic trajectory of the system. 

In this paper we show that it is possible to stabilize 
a periodic trajectory that preserves the qualitative fea- 
tures of a chaotic trajectory by injection of a periodic 
waveform. Our scheme is different from previous schemes 
of stabilization of a periodic trajectory in dynamical sys- 
tems 0) B 01 hi a way that in these previous examples the 
period of the stabilized trajectory is typically compara- 
ble to the time scales of the system dynamics, making the 
stabilized trajectory similar to a limit cycle more than to 
a chaotic trajectory. Here, on the other hand, the period 
of the stabilized trajectory is thousands of times longer 
than the typical time constants of the system dynam- 
ics and can be in principle arbitrarily long. This makes 
the stabilized periodic trajectory almost indistinguish- 
able from a truly chaotic trajectory of the system, and 
rigorously indistinguishable if the period of the control 
signal tends to infinity. We conjecture that the ensem- 
ble of periodic trajectories that may be stabilized has a 
continuous transition to the entire ensemble of chaotic 
trajectories as the period of the control waveform tends 
to infinity. Our specific example refers to the case of a 
semiconductor laser driven to chaos by optical feedback 
from a short optical cavity [5] . 

The dynamics of a semiconductor laser is described by 
the electric field of complex envelope E{t), normalized 
such that is the intra-cavity photon number, and 

by the carrier number N{t). These quantities are gov- 
erned by a modified set of Lang-Kobayashi (L-K) equa- 



tions fi'], which we generalize including an extra term 
accounting for the injection of an external driving field. 
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Here, No is the carrier number at transparency, G{t) is 
the saturated gain, g is the differential gain and s is the 
gain suppression coefficient, a is the line- width enhance- 
ment factor, Tp is the photon lifetime, Tc is the carrier 
lifetime, Tin is the round-trip time of the laser, rj is the 
fraction of the optical power fed-back to the laser cavity 
to drive the laser to chaos. At is the external cavity delay, 
ujo is the optical carrier angular frequency, / is the pump 
current and e is the electron charge. The term E(.xt{t) is 
the field amplitude of the injected external driving field 
after passing through the laser facet. Equations ([T]) and 
^ have been extensively used to study the laser dy- 
namics and, when compared with experimental results, 
they have passed the accuracy test in the totality of cases 
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We use first for the forcing field i?oxt(i) a replica of 
a portion of length T of the chaotic intra-cavity field of 
the free-running laser, attenuated by the factor ^/rf (rj' is 
the corresponding power fraction) , periodically repeated 
with period T. We find that, after an initial transient of 
the order of a few T, the intra-cavity field locks into a pe- 
riodic trajectory of period T. Results show stringent sim- 
ilarities with mode-locking. The periodic coherent injec- 
tion of a delayed waveform of the same laser acts like the 
coupling with an ideal external cavity of mode spacing 
1/T, which is a well known technique to achieve mode- 
locking with the period of the external cavity round-trip 
time T [8]. In this regime, the onset of a periodic evolu- 
tion may be seen as a form of chaos self-synchronization 
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FIG. 1: Intensity \E{t)f (dashed green curve) and \E{t-T)f 
(solid red curve), normalized to unit average value, of the 
electric field at time t and time t — T. Fig. [1^ refers to the 
free-running case. Fig. [TJ) to the locked case, where the two 
curves overlap on the scale of the plot. In Fig. [1]d, the blue 
dotted line represents the normalized intensity of the driving 
field at both time t and t — T. 




FIG. 2: Intensity \E{t)\'^ (dashed green curve) and \E{t — 
r)p (solid red curve), normalized to unit average value, of 
the electric field at time t and time t ~T, together with the 
normalized intensity of the driving field (dotted blue curve), 
for a) Ts = 25 ps, (b) Ta = 75 ps, and c) Ts = 150 ps sampling 
time. Circles denote the sampling points. 



[9|, or as a form of chaos control 0. In Fig. [T}i, we 
show by a solid red line the intensity of the field of the 
free-running laser, i.e. with £'ext(i) = 0, at time t — T, 
and with a dashed green line the same quantity at time 
t, normalized to unit average intensity. In Fig. [TJs, we 
report with a solid red line the intensity of the field at 
time t — T, with a dashed green line the same quantity 
at time t, and with a dotted blue line the intensity of 
the injected driving field i?ext(i), also normalized to unit 
average intensity. It is apparent that the laser intensity 
locks to the drive and into a periodic trajectory. We have 
verified that the full intra-cavity electric field, including 
the phase, and the carrier number are truly periodic, and 
that the periodic trajectory (once the laser parameters 
and the driving field are fixed) is a global attractor of 
the system, i.e. the system responds to a perturbation 
of electric field and carrier density returning to the same 
periodic trajectory. 

The values of the parameters of the numerical example 
are ^ Nq = 1.5 x 10^ 5 = 6 x lO-'^ns"!, s = 5 x 10"^ 
a = 4:, Tp — 5 ps, Tc — 0.5 ns. Tin = 8 ps, = 0.025, At — 
3 ns, uja = 2ttc/X with c the speed of light and A = 750 
nm, / = 1.5/th, with Ith = 60 mA. We choose 77' — 0.16, 
corresponding to a 16 % power coupling coefficient of the 
driving field. The period is T = 100 ns. 

In practice, the scheme described above may be re- 
alized in a master-slave configuration, by recording the 
complex chaotic waveform emitted by the free-running 
laser in a time interval T, and injecting in the same laser 
the output of a continuous-wave laser modulated in in- 
tensity and phase to periodically reproduce the recorded 
waveform. In real word, however, phase and amplitude 
modulation may only be applied within a finite band- 



width. For this reason, we decided to investigate the ef- 
fect on the locking quality of the modulation bandwidth. 
For this purpose, we applied to the slave a repeated copy 
of the free-running chaotic intra-cavity field sampled with 
a variable sampling time Ts- The sampled points were 
connected by an interpolation of the sampled intensity 
and phase using a MatLab interpolation routine based 
on spline (i.e. cubic) functions. The results are shown 
in Fig. [21 where we again compare the intensity of the 
field at time t — T with the intensity at time t and with 
the intensity, normalized as in Fig. (TJ of the forcing field 
E^^t{t), for T, = 25 ps, T, = 75 ps, and = 150 ps 
sampling time. 

We notice now that the laser locks into a periodic tra- 
jectory even when the sampling rate is lower than the 
bandwidth of the chaotic input, approximately 15 GHz in 
our numerical example. This is surprising because sam- 
pling the free-running chaotic waveform with a sampling 
rate significantly lower than its bandwidth is equivalent 
to selecting almost uncorrelated samples of the recorded 
chaotic trajectory. Stimulated by this observation, we 
decided to investigate a different configuration to achieve 
locking. A random phase modulation with uniform dis- 
tribution from to 2tt is applied with 1/Tm rate to a 
continuous-wave master, which is then injected into the 
slave. The intensity of the master laser is not modulated. 
We set the forcing field Ecxt{t) such that its intensity is 
equal to a fraction rj' of the average field intensity of the 
free-running slave. In Fig. [3] we show the results for 
T,„ ~ 25 ps, Tm ~ 75 ps, = 150 ps modulation time. 
For T,n = 25 ps the laser does not fully lock, although 
some periodic regularities show up. Full and robust lock- 
ing is instead achieved with Tm = 50 ps and larger. We 
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FIG. 3: Intensity \E{t)f (dasiied green curve) and \E{t-T)f 
(solid red curve), normalized to unit average value, of the 
electric field at time t and time t — T, together with the (con- 
stant) normalized intensity of the phase modulated driving 
field (dotted blue curve), for a) Tm — 25 ps, b) Tm = 75 ps, 
and c) Tm = 150 ps period of the random phase modulation. 



have verified in this case as well that, when the laser 
locks, the periodic trajectory is still a global attractor 
of the system. Although the curves of Figs. [2] and [3] 
show qualitative differences, simple inspection on a scale 
shorter than the period would identify both as plots of 
genuine chaotic time series. 

To obtain a quantitative characterization and a deeper 
understanding of the above observations, and to analyze 
the differences between sampled and random injection, 
we decided to geometrically characterize the attractor of 
the free-running laser by a numerical estimate of its Haus- 
dorff dimension D 1] . Our estimate is obtained using an 
A^-dimensional, with iV = 3 + 2At/St, approximation 
of the infinite-dimensional phase space, consisting of the 
value of N(t) and of samples of the real and imaginary 
parts of the electric field taken with St = 10 ps sampling 
time in the interval t — At < t' < t. In this reduced phase 
space, we follow the evolution of the system in a long time 
interval, which we chose as the time T, the period of the 
controlled regime. During a simulation run of T dura- 
tion, we periodically register with period t^ = 4:St a, full 
A^-dimensional picture of the dynamical system, ending 
up with a set of points in the A^-dimensional phase space 
(2500 points in our numerical example). We then nu- 
merically evaluate the correlation integral C(r), defined 
as the average number of A^-dimensional pairs with Eu- 
clidean distance less than r, averaged over all pairs, and 
we compute the exponent v characterizing the small r 
asymptotic power-law behavior of C(r) ~ r'^ . For the 
free-running laser, the chaotic attractor has Hausdorff 
dimension D c: 
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For a laser locked to a stable periodic trajectory the 
attractor is rigorously one-dimensional, but the asymp- 



FIG. 4: Exponent f for sampled injection (dots) vs. the sam- 
pling time Ts {Ts = 1 ps, ideal injection), and for random 
injection (squares) vs. the modulation period Tm- The value 
of 1/ of the free running laser is indicated with a star at Ts — 0. 
Dot-dashed line: number of estimated degrees of freedom of 
the input drive vs. the modulation period Tn for the case of 
random injection. In the inset, the curves of logj^Q [C'('')l vs. 
logjg (r) for the free-running (solid black) and the locked case 
with ideal injection (solid red). 



totic exponent v obtained using a time series of length T 
has a meaning by itself. If, for a sufficiently long period 
T, the exponent is close to the value of the standalone 
laser, we may reasonably conclude that the periodic tra- 
jectory covers a large part of the chaotic attractor of the 
free-running laser, and it would be hardly distinguishable 
from a genuine chaotic trajectory on a time scale of the 
order of T or shorter. In Fig. 2] we plot by red dots the 
values of ly vs. the sampling time T,, showing with a 
black star at = the exponent i/ of the free-running 
laser and with a red dot at = 1 ps that of the locked 
laser with ideal injection (i.e., sampling time equal to the 
numeric integration step, 1 ps). In the inset, we also plot 
the curves of log^g [C(r)] vs. logig(r) for these two cases, 
showing no difference between the two. This confirms 
that it would be impossible, over a time scale of the or- 
der of T, to distinguish the periodic trajectory from the 
truly chaotic one. These results are not surprising. In- 
deed, the laser locks with high fidelity to a drive obtained 
from an arbitrary segment of the laser electric field peri- 
odically injected into it. Consequently, the ensemble of 
all possible periodic trajectories will be arbitrarily close, 
for a sufficiently long period T, to a coverage of the en- 
tire periodic attractor. It might therefore be safely stated 
that the manifold of the ensemble of all the possible pe- 
riodic trajectories of the locked laser is the attractor of 
the free-running laser. This result holds whenever the 
control field is over-sampled, i.e. it is sampled with a 
sampling rate higher than its bandwidth, that is, in our 
case, for Tg < 75 ps. In this regime, Fig[5]shows that the 
exponent v is always equal to the free-running value. 
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For Ts > 75 ps, the free-running laser trajectory used 
as a drive is under-sampled. This means that the sampled 
input has less information than the free-running chaotic 
trajectory. To give an example, going from Tg — 75 ps 
to Ts — 150 ps the number of independent samples is 
halved, hence the input would reduce by half the num- 
ber of degrees of freedom. Figure |4] shows however that 
in this case the reduction by half of the number of de- 
grees of freedom of the input produces a reduction of the 
exponent v by only about 20 %. This may be an indica- 
tion of the property of the laser dynamics to recover part 
of the lost information of the drive, which is a natural 
trajectory of the system. 

Figure [4] reports also the values of i> for a random- 
phase modulated forcing field, plotted vs. with blue 
squares. With random injection, the random drive does 
not belong by any means to the free running chaotic 
attractor. In our A^-dimensional approximation of the 
phase space, the drive may be represented by points de- 
scribing in their temporal evolution a curve belonging to 
a manifold of dimension d < N. Indeed, by the Nyquist- 
Shannon sampling theorem with time and frequency in- 
terchanged, a time-limited complex signal in an interval 
At, which is a point in our phase space, is fully charac- 
terized by its spectral components (modulus and phase) 
spaced by AB = 1/ At. Although the number of spec- 
tral components to be included is virtually infinite, if the 
bandwidth of the drive is about 1/Tm, the number of 
effective degrees of freedom of the drive is d ~ At/Tm 
(half, because of the constant intensity, of the degrees 
of freedom of an intensity and phase modulated signal) . 
In the case of random phase modulation, the drive be- 
longs to a manifold of dimensionality d embedded into the 
phase space, which plays the role of the chaotic attractor 
of the sampled case. The ensemble of all possible laser 
outputs locked to the drive will belong to a manifold of 
dimensionality not larger than the dimension of the drive 
itself, because the laser, when locked, cannot create infor- 
mation. In Fig[4l we compare the exponent v estimated 
for the output with d* — 1.2 At/T„i, represented with 
a dot-dashed line, where the (empirically chosen) 20 % 
increment of d accounts for the slightly larger bandwidth 
of the input if compared to a rigorously band-limited sig- 
nal of bandwidth 1/Tm. The comparison shows that the 
d* ^ V hence that the laser output trajectory preserves 
all the degrees of freedom of the input. Notice that the 
value of V corresponding to Tm = 25 ps is below the 
l/r„i interpolation curve. This is because, in this case, 
the laser does not fully lock to the input, hence part of 
the degrees of freedom of the input are not reproduced 
by the output. 

With random injection, the manifold of the random 
drive docs not belong to the chaotic attractor of the free- 
running laser, hence the laser dynamics locked to the 
drive reproduce, at most, all the degrees of freedom of the 
drive itself. Consequently, the dependence of the expo- 



nent V on the bandwidth, hence on the number of degrees 
of freedom, of the drive is stronger here than with sam- 
pled injection. To give an example, reducing the degrees 
of freedom of the drive by half reflects into a proportional 
decrease of the value of exponent v of the output. 

To conclude, we have shown that the injection of a 
periodic waveform can lock a chaotic system into a peri- 
odic trajectory belonging to its attractor. The resulting 
process is periodic with a local chaotic signature, with a 
dimensionality that is different from the standalone pro- 
cess, depending on the specific injection scheme. Our 
analysis is specialized to a system obeying the modified 
L-K equations, appropriate, for example, to the descrip- 
tion of an injected semiconductor laser driven to chaos 
by optical feedback from a short external cavity. Here, 
the injection of a periodically replicated sampled copy of 
the laser output and a continuous- wave randomly-phase- 
modulated optical field lock the laser into a process that 
is locally equivalent to a chaotic dynamic. A careful anal- 
ysis of the locked output reveals that with sampled injec- 
tion the output waveform laser covers the same attractor 
of the standalone laser, whereas with random injection 
it belongs to a manifold of dimensionality equal to the 
number of degrees of freedom of the drive. 
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